Department  of 
Electrical  and  Computer 

Engineering 

Bainesville,  Florida  32611 


UNIVERSITY  OF 

FLORIDA 


20071 1 1 6225 


MIMO  Radar 

-  Diversity  Means  Superiority 

Jian  Li 

Department  of  Electrical  and  Computer  Engineer  inn 
P.O.Box  116130 

University  of  Florida,  Gainesville,  FL  32611 
Phone:  (352)  392-2642;  Fax:  (352)  392-0044;  Email:  li@dsp.ufi.edu 

Annual  report  for  the  Office  of  Naval  Research 
Grant  No.  N00014-07- 1-0293 
November  2006  -  October  2007 


Best  Available  Copy 


REPORT  DOCUMENTATION  PAGE 

Form  Approved 

OMB  No.  0704-01 88 

Th*  public  r apart  mg  burden  for  this  collection  of  information  is  aslirnatad  to  avarage  1  hour  par  rasponse.  including  tha  time  tor  reviewing  instruction*.  starching  axuting  dais  sources. 

flathiring  and  maintaining  tha  data  n«*dadr  and  completing  and  reviewing  th*  collection  of  information.  Sand  comments  regarding  Hus  burdan  estimate  or  any  other  aspect  of  this  collection  of 
information,  meludmg  suggestions  for  reducing  the  burden.  to  Department  of  Defanse.  Washington  Haadqu after S  Service*.  Dirac  1  orate  tor  Information  Operations  and  Reports  (Q7CH  OlflB' 
121$  Jefferson  Davis  Highway.  Suita  1204.  Arlington,  VA  22202-4302.  Respondents  should  be  aware  That  ho twilhst ending  any  other  provision  of  law.  no  parson  shall  ba  sublet  to  any 
panel ty  fnr  failing  io  comply  with  a  collection  of  information  if  it  does  not  display  a  currently  valid  OMB  control  number 

PLEASE  DO  NOT  RETURN  YOUR  FORM  TO  THE  ABOVE  ADDRESS, 

1.  REPORT  DATE  (DD-MM-YYYY1 
25-10-2007 

2  REPORT  TYPE 

Annual  Report 

3.  DATES  COVERED  IFrom  Toi 

Nov.  2(ll)(i  -  (ki  20117 

4.  TITLE  AND  SUBTITLE 

Multi-Input  Multi-Output  (MlMO)  Radar  -  Diversity  Means  Superiority 

5a,  COr\ 

(TRACT  NUMBER 

N000 14-07-1 -0293 

5b.  GRANT  NUMBER 

N0UU 14-07- 1-1)293 

5c.  PROGRAM  ELEMENT  NUMBER 

6,  AUTHORS 

Li,  Jian 

5d.  PROJECT  NUMBER 

00064512 

5e.  TASK  NUMBER 

51,  WORK  UNIT  NUMBER 

7,  PERFORMING  ORGANIZATION  NAMEIS)  AND  ADDRESS(ES) 

Department  of  Electrical  and  Computer  Engineering 

P.  O.  Box  1 16130,  University  of  Florida 

Gainesville,  FL  32611 

8  PERFORMING  ORGANfZATION 

REPORT  NUMBER 

N  A 

9,  SPONSORING/MONITORING  AGENCY  NAMEIS1  AND  ADDRESSEES) 

Office  of  Naval  Research 

S75  North  Randolph  Street 

Arlington,  VA  22203-1995 

10,  SPONSOR/MONITOR'S  ACRONYMIS) 

ONR 

11  SPONSOR/MONITOR  s  REPORT 

NUMBERIS) 

\  \ 

12,  DISTRIBUTION /AVAILABILITY  STATEMENT 


Approved  for  public  release;  distribution  unlimited 

13,  SUPPLEMENTARY  NOTES 

None 

14.  ABSTRACT 

We  consider  a  multiple-input  multiple-output  (MlMO)  radar  system  where  both  the  transmitter  and  receiver  have  multiple  wcl 
separated  subarrays  with  each  subarray  containing  dosely*spaecd  antennas.  Because  of  this  general  antenna  eonflgunumn.  both  t In¬ 
coherent  processing  gain  and  the  spatial  diversity  gain  can  be  simultaneously  achieved  We  compare  several  spatial  speuul 
estimators,  including  Capon  and  APES,  for  target  detection  and  parameter  estimation.  We  introduce  a  generalized  iiMtinhM  i  n  . 
test  (GLRT)  and  a  conditional  generalized  likelihood  ratio  test  (cGLRT)  for  the  general  antenna  amgu ration  Based  on  K  d  K  I  and 
cGLRT,  we  then  propose  an  iterative  GLRT  (iGLRT)  procedure  for  target  detection  and  parameter  esimiauon  \  ia  several 
numerical  examples,  we  show  that  iGLRT  can  provide  excellent  detection  and  estimation  performance  at  a  low  compuianms.il  on 


15,  SUBJECT  TERMS 

Adaptive  arrays,  MlMO  Radar,  detection,  localization,  parameter  estimation,  growth-curve  (GC)  model,  block-diagonal  growth- 
curve  (BDGC)  model,  generalized  likelihood  ratio  test  (GLRT). 


16.  SECURITY  CLASSIFICATION  OF; 

17.  LIMITATION  OF 

ABSTRACT 

18.  NUMBER  l 

19a,  NAME  OF  RESPONSIBLE  PERSON 

a.  REPORT 

b,  ABSTRACT 

c.  THIS  PAGE 

OF 

PAGES 

Jian  Li 

U 

U 

u 

U 

39 

19b.  TELEPHONE  NUMBER  Unclurte  area  rortt" 

(.352)  392-2642 

Standard  Form  298  iRifv  R  ‘<8 

^flsr.i  'btfl  n  AMS’  Stti  '  i  ^  - 


TABLE  OF  CONTENTS 


MW 

LIST  OF  FIGURES . . .  Hi 

ABSTRACT  . .  H 


CHAPTER 

1  INTRODUCTION .  I 

2  SIGNAL  MODEL .  1 

3  SPATIAL  SPECTRAL  ESTIMATORS .  (i 

3.1  Capon .  U 

3.2  APES .  7 

-1  GENERALIZED  LIKELIHOOD  RATIO  TEST .  ') 

LI  Generalized  Likelihood  Ratio  Test.  [GLUT] .  < 

1.2  Conditional  Genera  lined  Likelihood  Ralio  I  >  >t  it.  I  I 


•1.3  Iterative  Generalized  Likelihood  Ratio  Test  (iCI.lil  i *>■ 

5  NUMERICAL  EXAMPLES .  21) 

5.1  Cramer-Rao  bound .  ill 

5.2  Target  Detection  and  Localization .  22 

6  CONCLUSIONS .  25) 

A  CRAMER-RAO  BOUND  .  3(1 

REFERENCES . 33 


Best  Available  Copy 


LIST  OF  FIGURES 


Figure  page 

5  1  Cumulative  density  functions  of  the  Cmtner-Rao  bounds  for  (a)  0,  .mi] 

(b)  ftt.  21 

5  2  Outage  CRB  versus  SNR.  fa)  CUBo.ih  for  0,,  (b)  C'RBu.ni  t'«»i  ">  » 

CRB,,. i  for  Q\,  mid  (d)  CRB0,,  for  ft, .  21 

5  3  Spatial  spectra,  and  GLR  and  cGLR  Pseudo-Spectra  .when  0\  It) 

02  —  -20p,  mid  -  0".  (a)  LS.  (b)  Capon,  (cj  APES.  (di  t . I . I M  .md 
(e)  iGLRT . . . . 

5  1  Spatial  spectra,  and  Gldi  and  cGLR  Psciido-Spe<  i  m.  when  » 

02  ~  —  P,  and  ftv=t)°.  (a)  Capon,  (b)  APES.  (<■)  GUI  1  and  ,i  a  .  I ;  I  2o 

5  5  GLR  mid  cGLR  Pseudo-Spec  tru  obtained  in  Steps  I  mid  II  <it  <;|.RJ 
when  0\  =  — 4QC ,  On  -  -4°,  and  0.-,  -  ON  fa)  f>(0).  flu  />{(/[()  < 

p{0\0\Jhi).  and  (d)  p(0|01;  ft,.  ftR .  27 

5  ti  Cumulat  ive  density  functions  of  the  CRBs  and  MSEs  for  (a)  0\  and  (b)  ftt  27 

5  7  Outage  CRBiu  and  MSEn.i  versus  SNR  for  (a)  0\  and  (b)  ftt .  28 

5  8  Outage  CRB(,,i  and  MSEn.i  versus  number  of  data  sampli*  I'm  .0  0  ,.i  d 
(b)  ft, . 


iii 


Abstract 

We  consider  a  multiple-input  multiple-output  (MIMO)  radar  system  where  Imili  i  h< 
transmitter  and  receiver  have  multiple  well-separated  subarrays  with  each  MtUm.n 
containing  closely-spaced  antennas.  Because  of  this  general  mitemin  <  i m h<u m . l 1 1 >  > r > 
both  the  coherent  processing  gain  and  the  spatial  diversity  gain  can  he  sinuili  aue- 
ously  achieved.  We  compare  several  spatial  spectral  estimators,  including  Capon 
and  APES,  for  target  detection  arid  parameter  estimation.  We  introduce  a  general¬ 
ized  likelihood  ratio  test  (GLRT)  and  a  conditional  generalized  liki  lilmud  i.nm  i<  -i 
(cGLRT)  for  the  general  antenna  configuration.  Based  on  GLRT  and  cGLR'l .  wi 
then  propose  an  iterative  GLRT  (iGLRT)  procedure  for  target  detection  mid  param¬ 
eter  estimation.  Via  several  numerical  examples,  we  show  that  iGLRT  can  provide 
excellent  detection  and  estimation  performance  at  a  low  computational  cost. 
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CHAPTER  1 
INTRODUCTION 

A  multiple-input  multiple-output  (MIMO)  radar  uses  multiple  antemms  in  si¬ 
multaneously  transmit  several  linearly  independent  waveforms.  It  also  uses  multiple 
antennas  to  receive  the  reflected  signals.  It  has  been  shown  that  by  exploiting  Mils 
waveform  diversity,  MIMO  radar  can  overcome  performance  deg rad at  ions  <  a  used  i  In 
radar  cross-section  (RCS)  fluctuations  [i]  -  [4],  achieve  flexible  spatial  transmit  beam- 
pattern  design  [fj]  [6],  provide  high-resolution  spatial  spectral  esiimnt.es  .uni 

significantly  improve  tiie  parameter  identifiability  s] 

The  statistical  MIMO  radar,  studied  in  [l]  -  [I],  aims  ni  resisting  tin  -i  miill.i 
tion”  effect  encountered  in  radar  systems.  It  is  well-known  that  the  IK'S  u!  .i  i,n  .  i 
which  represents  the  amount  of  energy  reflected  from  the  target  toward  the  m-rivei 
changes  rapidly  as  a  function  of  the  target  aspect  [19],  and  the  locations  of  t  he  trans¬ 
mitting  and  receiving  antennas.  The  target  scintillation  causes  severe  degradations  m 
the  target  detection  and  estimation  performance  of  the  radar.  By  spacing  the  trans¬ 
mit  antennas,  which  transmit  linearly  independent  signals,  far  away  from  each  other, 
a  spatial  diversity  gain  can  be  obtained  as  in  the  MIMO  wireless  cotunmnicatums  to 
this  scintillation  effect  [1]  -  [1]. 

Flexible  transmit  beampattern  designs  are  investigated  in  [a]  [(>],  Different  h  um 
the  statistical  MIMO  radar',  the  transmitting  antennas  are  closely  spaced.  The  au¬ 
thors  in  [0]  and  [(>]  show  that  the  waveforms  transmitted  via  closely  spaced  antennas 
can  be  optimized  to  obtain  several  transmit  beampattern  designs  witli  siijnrtoi  pi-i 
formance.  For  example,  the  covariance  matrix  of  the  waveforms  can  be  optimized 
to  maximize  the  power  around  the  locations  of  interest  and  also  to  minimize  the 
cross-correlation  of  the  signals  reflected  back  to  the  radar  by  these  targets,  t  herein 
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significantly  improving  the  performance  of  the  adaptive  MIMO  radar  techniques.  Due 
to  the  significantly  larger  number  of  degrees-of- freedom  of  a  MIMO  system,  impmvrd 
transmit  beampatterns  can  be  achieved  with  a  MIMO  radar  than  with  its  phased- 
array  counterpart. 

In  [!)],  a  MIMO  radar  technique  is  suggested  to  improve  the  radar  resolution 
The  idea  is  to  transmit  N  (Ar  >  1}  orthogonal  coded  waveforms  by  N  antennas  and 
to  receive  the  reflected  signals  by  M  (M  >  1)  antennas.  At  each  receiving  anteiimi 
output,  the  signal  is  matched-filtered  using  each  of  the  transmitted  wavelorms  in 
obtain  NM  channels,  where  the  data-adaptive  Capon  beamformer  [2d]  is  applied,  Ii 
is  proved  in  [Si]  that  the  beampattern  of  the  proposed  MIMO  radar  is  obtained  la¬ 
the  multiplication  of  the  transmitting  and  receiving  beampatterns,  hence  it  has  high 
resolution.  However,  [d]  considers  only  single-target  scenarios. 

A  MIMO  radar  scheme  is  considered  in  [lfi]  and  [IT]  that  can  deal  with  the 
presence  of  multiple  targets.  Similar  to  some  of  the  MIMO  radar  approaches  [7]  -  'IT), 
linearly  independent  waveforms  are  transmitted  simultaneously  via  multiple  antennas 
Due  to  the  different  phase  shifts  associated  with  different  propagation  paths  fmm 
transmitting  antennas  to  targets,  these  independent  waveforms  are  linearly  combined 
at  the  targets  with  different  phase  factors.  As  a  result,  the  signal  waveforms  reflected 
from  different  targets  are  linearly  independent,  allowing  the  direct  application  of  mam- 
adaptive  techniques  to  achieve  higher  resolution  and  interference  rejection  eapabiliiv 
Several  adaptive  nonparametric  algorithms,  some  of  which  also  model  steering  vet  on 
errors,  are  presented  in  [Hi]  and  [17], 

The  MIMO  radars  discussed  above  can  be  grouped  into  two  classes  according  to 
their  antenna  configurations,  One  class  is  the  conventional  radar  array,  in  which  both 
the  transmitting  and  receiving  antennas  are  closely  spaced  for  coherent  transmission 
and  detection  [5]  -  [17],  The  class  other  is  the  diverse  antenna  configuration,  where 
the  antennas  are  separated  far  away  from  each  other  to  achieve  spatial  dio-i.'in  g.mi 


[]]  -  [1].  To  exploit,  the  benefits  of  both  schemes,  we  considei  a  geurial  . 1 1 ti <  1 1 1 1. 1 
configuration  in  this  report,  i.e.,  both  the  transmitting  and  receiving  antenna  nunv* 
consist  of  several  well-separated  subarrays  with  each  subarray  containing  cluseb 
spaced  antennas.  We  establish  the  growth  curve  models  [21]  -  [21]  and  devise  several 
estimators  for  the  proposed  MIMO  radar  system. 

The  remainder  of  this  report  is  organized  as  follows.  Chapter  .  presents  tin1 
MIMO  radar  signal  model.  In  Chapter  .'1,  we  discuss  two  adaptive  spatial  spectral 
estimators  including  Capon  [20]  and  APES  [25].  In  Chapter  1,  we  introduce  a  gen¬ 
eralized  likelihood  ratio  test  (C1LRT)  and  a  conditional  generalized  likelihood  ratio 
test  (cGLRT),  and  then  propose  an  iterative  GLRT  (iGLRT)  procedure  for  target 
detection  and  parameter  estimation.  Numerical  examples  are  provided  in  Chaptci 


CHAPTER  2 
SIGNAL  MODEL 

Consider  a  narrow-hand  MIMO  radar  system  with  N  and  M  subamiys  fm  i  mns- 
mitting  and  receiving,  respectively.  The  nth  transmit  and  mth  receive  suharrays 
have,  respectively,  Nn  and  Mm  closely-spaced  antennas,  n  =  1,  2,  .  Y.  and 

m  =  1,  2,  ■■■  ,  M.  We  assume  that  the  subarrays  are  sufficient l\  sej  wanted  .nn! 
hence,  for  each  target,  its  radar  cross-sections  (RCSJ  for  .different  t  ninsmu  ,iml  n  <  n  . 
subarray  pairs  are  statistically  independent  of  each  other,  Let  vr, (9 )  and  a,,, ( tl  i  la 
the  steering  vectors  of  the  nth  transmitting  subarray  and  the  mth  receiving  sulmmiv 
respectively,  where  0  denotes  the  target  location  parameter,  for  example  its  angular 
location.  Let  the  rows  of «!»,,  be  the  waveforms  transmitted  from  the  antennas  of  the 
nth  transmit  subarray.  We  assume  that  the  arrival  time  is  known.  Then,  the  signal 
received  by  the  mth  subarray  due  to  the  reflection  of  flit'  target  at  0  can  la  w  linm 
as 

N 

X"’  =  Vj (0)$„  +  Z,n,  III =  1.  ■  ■  ,  ,  J  |  i 

n=l 

where  0mnto  is  the  complex  amplitude  proportional  to  the  HCS  tot  the  [m.  //lilt 
receive  and  transmit  subarray  pair  arid  for  the  target  at  the  location  0  The  matrix 
ZT„  denotes  a  residual  term  for  the  unmodelled  noise,  e.g.,  interferences  Imm  t.irget- 
other  than  0  and  at  other  range  bins,  and  intentional  or  iminteutiiiiijil  |,imnm.  ■ 
notational  simplicity,  we  will  not  show  explicitly  the  dependence  oi  Z,„  mi  w 
Let 


X=[X[  ■••Xr/6C^, 

(2  2) 

A(0)=Diag[a,(0),  ,  a* (*)]  e  C"xA>, 

(2  3) 

V(0)  =  Diag[V|(0),  ■■■  ,vt,{9)\€CN**. 

(2  1) 

4 


and 


O  =  [0'f  $^,]r  E  C,v*'\ 

where  M  =  M\  +  •••  +  M ^  and  N  —  N\  +  +  N#  are  the  total  lmnibers  ul 

receive  and  transmit  antennas,  respectively,  L  is  the  number  of  data  samples  of  the 
transmitted  waveforms,  (-)r  denotes  the  transpose  operator,  and  Diaglai,  •  ■  •  aw  I  i> 
a  block-diagonal  matrix  with  at,-*  -  , being  its  diagonal  sitlmuUritvs  11  n  n 
can  be  rewritten  in  the  growth-curve  (GC)  model  [‘J I ]  -  [22]: 

X  =  A(0)BeS(0)  +  Z,  (2  b) 

where  the  (m,  n)th  element  of  the  M  x  N  matrix  B#  is  Z  is  defined  similarly  to 

X  in  (2  2),  and  the  rows  of  S{9)  are  the  reflected  waveforms  by  the  target  at  location 
0,  i.e., 

S(0)  =  Vr(fl)<i>.  j: 

Note  that  when  N  —  M  —  1,  the  signal  model  in  (2  ii)  ret  burs  to  the  MIMt  ) 
radar  model  in  [ J (>]  -  [17],  whereas  when  N  =  N  and  M  =  M  it  rctlttco  to  tin 
diversity  data  model  in  [1]  -  [.*.{].  Based  on  this  data  model.  We  below  derive  two 
classes  of  nonparametric  methods,  i.e.,  spatial  spectral  estimation  and  generalised 
likelihood  ratio  test  (GLRT),  for  target  detection  and  localization. 


CHAPTER  3 

SPATIAL  SPECTRAL  ESTIMATORS 

We  discuss  two  spatial  spectral  estimators  for  the  proposed  MIMO  radm  system 
We  use  these  methods  to  estimate  the  complex  amplitudes  in  B„  for  cadi  H  <>!  mi  not 
from  the  observed  data  matrix  X.  The  Frobeuius  norm  of  the  estimated  Q„  forms 
a  spatial  spectrum  in  the  ID  case  or  a  radar  image  in  the  2D  case.  We  ran  t  hen 
estimate  the  number  of  targets  and  their  locations  by  searching  foi  the  peaks  m  i|« 
estimated  spectrum  (or  image). 

A  simple  way  to  estimate  Bo  in  (2  fi)  is  via  the  Least-Squares  (LSI  method 

B  1.3,0  =  [A"(0)A(e)}-iA(0)XSH(9)[S(B)S"(0)}-'.  ,3  I 

where  (  )t!  denotes  the  conjugate  transpose.  However,  as  any  other  data-independent 
beamforming-type  method,  the  LS  method  suffers  from  high-sidelobes  and  low  reso¬ 
lution.  In  the  presence  of  strong  interference  and  jamming,  the  method  completely 
fails  to  work.  Hence,  we  discuss  two  robust  adaptive  spatial  spectral  estimation 
approaches  that  offer  higher  resolution  and  interference  suppression  capabilities. 

3.1  Capon 

The  Capon  estimator  for  By  in  (2  (j)  consists  of  two  main  steps  [20],  [2fi],  [22]. 
The  first  step  is  a  generalized  Capon  beamforming  step,  The  second  step  is  an  LS 
estimation  step,  which  involves  a  matched  filter  to  the  known  waveim-m  S1  c 
The  generalized  Capon  beanaformer  can  be  formulated  as  follows 

rnintr(WHRW)  subject  to  W11  A(ti)  =  I.  (32) 

where  W  6  C MxM  is  the  weighting  matrix  used  to  achieve  noise,  interference  and 
jamming  suppression  while  keeping  the  desired  signal  undistorted,  tr(-)  denotes  the 


0 


trace  of  a  matrix,  and 


R  =  jXXH 

Lj 


(3  3) 


is  the  sample  covariance  matrix  with  L  being  the  number  of  data  samples. 

Solving  the  optimization  problem  in  (3  2},  we  have: 

Wc^^R-UCejlA^^R-'Ate)]-1.  (3  1) 

By  using  {2  (5)  and  (3  I),  the  output  of  the  Capon  bearnfornier  can  be  wiirteu 

iA/J(^)R-1A(^]’IA;/(^)R-1X  =  B„S(6/)  +  [Atf(fl)R"'A(fl)]-IA/,(fl)R  *2.  (3  5] 

By  applying  the  least-squares  (LS)  method  to  (3  5),  the  Capon  estimate  of  follows. 

Bcw  =  [A//(0)R.-1A(0)]-1A^(0)R.-1XS"C@)[S(0)Swr(0)J-1.  (3-6) 

3.2  APES 

The  generalized  APES  method  is  a  straightforward  extension  of  t  he  A PES  met  hud 
[25]  [27],  which  can  be  formulated  as: 

min  ||  \\'nX  -  B«S(0)  ||J  subject  to  WwA(0)  =  1  '17' 

where  ||  •  ||  denotes  the  Frobenius  norm,  and  W  is  the  weighting  matrix.  Minimizing 
the  cost  function  in  (3  7)  with  respect  to  Bfl  yields: 

BAPES,e  =  W"XS"(0)[S(0)S"(tf)l-'.  (3  S) 

Then,  the  optimization  problem  reduces  to 

mintrfW^QW)  subject  to  W/yA($)  =  I,  (3  9) 

with 

Q  =  R-  jxS"(0)[S(0)S"(fl)]"‘S(tflX/V  ,:i  III 

For  notional  simplicity,  we  have  omitted  the  dependence  ol  Q  ou  0 
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Solving  the  optimization  problem  oi'  (3  !))  gives  the  generalized  APES  beam 
former  weighting  matrix: 

WAPES)9  =  Q_1A(0)[Aw(i9)QA{0)]-1,  (3-  11} 

Inserting  (3  11)  in  (3  8),  we  readily  get  the  APES  estimate  of  as: 

BAPB3,fl  =  [AH(6)Q-lA(0)]-lAH(e)Q-lXSH($)[S(Q)SH{e)}-1.  (3  12) 

Interestingly,  we  note  that  (3  12)  has  the  same  form  as  the  ML  estimate  in  [21] 
and  [22],  However,  the  APES  estimate  is  derived  based  on  the  beamforming  method, 
and,  unlike  the  ML  in  [21]  and  [22],  it  does  not  need  probability  density  function 
(pdf)  of  Z. 


CHAPTER  4 

GENERALIZED  LIKELIHOOD  RATIO  TEST 


Generalized  likelihood  ratio  test  (GLRT)  has  been  used  widely  foj  target  de¬ 
tection  and  localization.  We  derive  below  a  GLRT  and  a  conditional  genei  aliped 

likelihood  ratio  test  (cGLRT)  for  the  proposed  Ml  MO  rndni.  and  then  pi>i| . 

iterative  GLRT  (iGLRT)  procedure  for  improved  performance. 

4.1  Generalized  Likelihood  Ratio  Test  (GLRT) 

Throughout  this  chapter,  we  assume  that  the  columns  of  the  interference  and 
noise  term  Z  in  (2  G)  are  independently  and  identically  distributed  (i.i.d.)  circu¬ 
larly  symmetric  complex  Gaussian  random  vectors  with  mean  zero  and  an  unknown 
covariance  matrix  Q, 

Consider  the  following  hypothesis  test,  problem: 


J  Hn  :  X  -  Z 

H]  :  X  =  A(0)BflS(iP)  +  Z. 


i  I 


i.e.,  we  want  to  test  if  there  exists  a  target  at  location  8  or  not  Similarly  to  Nj  and 
[29],  we  define  a  generalized  likelihood  ratio  (GLR)  as  follows: 


pm  = 


}t 


I  2 


where  /(X]Hj)  (i  —  0,  1)  is  the  pdf  of  X  under  the  hypothesis  Ht.  From  (  1 ).  we 

note  that  the  value  of  the  GLR,  p{8),  lies  between  0  and  1.  If  there  is  a  target  at 
a  location  0  of  interest,  we  have  rnaxe  q  3>  maxq i.e.,  p  %  1: 

otherwise  p  ~  0, 
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Under  Hypothesis  Ho,  we  have 

/(xlHo)  =  exp(-  trlQ-^XX")},  (4  :i) 

wiiere  j  ■  |  denotes  the  determinant  of  a  matrix.  Maximizing  ( 1  3)  with  respect  to  Q 
yields: 

max/(X|Ho)  =  (7re)_iM|R|"/',  (I  li 

where  R  is  defined  in  (3  3). 

Similarly,  under  Hypothesis  H  ( ,  we  have 

f(X\Hl)  =  ^n^exp{-tt{Q-1[X- A($)B0S(e)][X- A(d)Bt,S(ff)^}} 


Maximizing  (  I  D)  with  respect  to  Q  yields: 


max/(X|Hj)  =  (ttc) 


~LM  ^[X  -  A(6/)BeS(0)][X  -  A(0)BoS(0)]" 


I  I  l-i 


Hence,  the  optimization  problem  in  the  denominator  of  (1  1)  reduces  to 


mm 

Bo 


-\X  -  A(0)BeS(0)][X  -  A(0)BoS(fl)]" 


(4  7) 


Following  [21]  and  [22]  and  dropping  the  dependence  of  A.  S  anti  B  ou  0  fuj 


notional  convenience,  we  have: 


j[X-  ABS][X  -  ABS]" 

X/ 

Y  [AB  -  XSw(SS")-*](SSh)[AB  -XS"(SS"rl]W  +  Q 
=  |Q||i+|q-*[AB  -XS"(SSh)"1](SS")[AB-XS"(SSh)"1]wQ  *  I 

=  |Q  1+  }(SS")s[AB  -  XS^SS")-1]"^-1  X 

J-i 

[AB-XS"(SS")-t](SS")4 

-  IQI  1  +  2{ssh)-4sxw[Q-'  -  Q-,A(A"Q-|A)-|A,'Q  ']XS"[SS" 

/  j 

y(SS")5[B  {AnQ~  1  A}- 1  AhQ_1XSw(SSh)  ']"(A"Q  'Al  * 

[B  -  (A//Q"1A)“lA"Q^‘XS"(SSw)“t](SS//)5 
>  |Q|  1+  y{SS")-5SXw  [Q-1  -^"1A(A"Q-,A)-IAffQ’I]XSiE^SS")t4 

i_ t 


9) 


where  Q  is  defined  in  (3  10).  To  get  {!  K),  we  have  used  the  fact  that  |[  +-  XY 
|I  +  YX|  [311],  and  the  equality  in  (1  0)  holds  when  B  equates  to  the  APES  estimate 
in  (3  12).  Note  that 


IQI 


I  +  -(SS")“5SX"  [Q-1  *  Q-IA(A"Q-1A)_I  AHQ-']XSW(SS,/)  s 

Zj 


=  IQI  |I  +  [Q-1  -  Q“J  AfA^Q-'A)-1  A^Q-^R  -  Q)| 
-  |R  —  A(AWQ_1  A)~*  AwQ-l(R  —  Q)| 

=  [R||I  -  A//(A//Q“lA)-1A“(Q"1  -  R  '}| 

=  |R||(a"q-1a)-1(a"r-1a}|1 


I  1  In 


From  (  I  (i),  (  I  9)  and  (  !  10),  it  follows  that 


max /(X|Hi)  =  (ffe)"t"|R[-t|A"(ejQ"1A(fl)|i|A"(0)R'lA(^|"t.  (1  11) 

t3n,Q 
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Substituting  (1  1)  and  (1  II)  into  (1  2)  yields: 


m  = 


We  remark  that  when  there  are  multiple  targets,  and  the  unmbei  ut  tuigeis 
(say  K)  are  known  a  priori,  the  GLRT  in  (  I  12)  can  be  extended  to  a  multivariate 
counterpart  by  considering  the  following  hypothesis  testing  problem: 


(4-13) 


As  a  parametric  method,  this  multivariate  GLRT  can  provide  better  target  detection 
and  parameter  estimation  performance  than  its  univariate  counterpart  I  lowevei  the 
multivariate  GLRT  is  computationally  intensive  because  it  needs  to  search  in  tin 
A' -dimensional  parameter  space  {0fc}f=1.  Moreover,  the  number  of  targets  is  ha  id  Is 
known  a  priori  in  practice. 

We  propose  below  an  iterative  GLRT  (iGLRT),  which  require  only  one-dimensional 
search  (like  the  univariate  GLRT),  but  provides  a  target,  detection  and  parameiei  es¬ 
timation  performance  close  to  the  multivariate  GLRT, 

4.2  Conditional  Generalized  Likelihood  Ratio  Test  (cGLRT) 

Before  we  describe  the  iGLRT  procedure,  we  first,  consider  the  following  hypoth¬ 
esis  testing  problem,  referred  to  as  the  conditional  generalized  likelihood  ratio  lest 
(cGLRT).  Suppose  that  we  know  that  there  are  p  targets  at  the  locations  { Bk  | . 
and  we  want  to  determine  if  there  are  any  additional  targets.  This  problem  can  be 
formulated  in  the  following  hypothesis  testing  problem: 


HP  :  X-S*_|A(0*)B^S(^)  “  Z 

HJ)+1  :  X  -  A(0)BflS(0)  +  XX,  Add  !  It..  SRd  )  Z 
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Note  that  both  the  equations  in  (1  1  I)  are  ill  the  form  of  the  block-diagonal  growth 
curve  (BDGC)  model  studied  in  [21],  For  convenience,  we  rewrite  (  I  II)  as: 

f 

Hj,  :  X  AjjBjjSjj  +  Z 

<  ...  (i  to) 

Hjn-j  .  X  Ap+iBp+iSp+i  +  Z, 

where 


B, 

=  DiagfB,^.  ■■■  . 

Bu  ) 

< 

K 

II 

'> 

> 

s 

i  i; 

sP 

=  [Sr(0i)  Sr 

i  i** 

Bp*i 

=  Diag(Bw,B^, 

■■■ > 

!  1  PJ 

Ap+i 

=  [A(6I)  A(^)  •• 

■  A(0P)], 

(  i  20) 

Sp+1 

=  [Sr(0)Sr(0,)  ■ 

■■■  sr(^)f. 

(  !  21) 

Similarly,  we  define  a  conditional  generalized  likelihood  ratio  (cCLR)  as  follows 


maxBr,.cj  /(X|HP) 

1  \ 
t 

maXBp+j,Q  /(X|Hp+[ ) 

/ 

(  I  221 


where  /(X|H,)  is  the  pdf  of  X  under  the  H,  hypothesis,  and  Q  is  the  covariance 
matrix  of  the  columns  of  Z. 

We  first  consider  the  optimization  problem  of  the  mmimitor  i:i<  M,i\i 

mizing  /{X|HJt)  with  respect  to  Q  yields: 


max /(X'Hp)  =  (7re)~ML 


T(X  -  A,,B(1Sft)(X  -  A,,BpSp)H 

JU 


(  1  23  j 


Hence,  the  optimization  problem  reduces  to 


min 

b„ 


— (X  -  ApBpSp)(X  —  ApBpSp)^ 


with  BfJ  -  Diag(B(j|,  *•■  ,  B^).  (4  24) 


1-1 


The  optimization  problem  in  (I  'J  i)  does  not  appear  to  admit  a  closed-form  solution 
because  Bp  is  a  block-diagonal  matrix.  Herein,  we  adopt  a  technique  used  in  j;nj  to 
approximate  a  closed-form  solution. 

Note  that 


where 


and 


|(X  -  A„B„S„)(X  -  AfJB,,S„)" 

^-(x  -  ApBpSpXn^,  +  n±  )(x  -  ApBpS„)" 

i-(xns;j  -  ApBpSp) (Xllg^  -  ApBpSp)"  +  Q„ 


f  (xn§ii  -  ApBpSp) (xn^  -  ApBpSp)" q;1  + 1 


(i  -h  ) 


l\,  =  S"(SpS")-Sp, 


=  i-ru 


(  I  2(\\ 


Q„  =  -xn|x", 


I  I  271 


with  (■)“  denoting  the  generalized  matrix  inverse. 

Consider  the  idempotent  matrices  rig  and  .  Assume  that  the  inmibn  nt 
data  samples  is  large  enough,  i.e.,  L  3>  Np.  Note  that  S,,  is  an  A  />  x  I.  maim 
Hence,  we  have: 

rank(IIg(i)  <  Np  and  rank  (11^  )  >  L  —  Np ,  (4  28) 


with  rank(-)  denoting  the  rank  of  a  matrix.  Then,  we  have 


Q„  =  O(i) 


(1  29) 
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and 

~(XU^r  -  ApBpSp)(XIIsfi  -  A,,B;jSp)r 

=  -^{X  -  ApBpSp)n§p(X  -  ApBpSp)7,  =  ^  (7)  ■  (  I  30) 

Therefore,  we  get: 

7(xn§p  -  ApBpSp) (xn§;i  -  ApBpsp)TQJ71  =  o  Q-)  «  1.  (4 -31) 

Let  {A,  }?li  be  the  eigenvalues  of  the  matrix  iu  (  I  11 ),  which  satisfy  e  hat  II  \  l 
Through  some  matrix  manipulations,  we  obtain: 


(XLIgf  —  A,3;,S?J)(Xn^  —  ApBpSfJ)/  QfJ  ]  4-  I 


M 


M 


=  JJ(i  +  y) « 1  +  53  a, 


i= 1 


1  7 


1+Ltr 

1  +  7  !l  vec(Qp  *Xngp)  -  vec(q; 7  ApBpSp)  ||2, 


(xn§p  -  j ApBpSp)(Xngp  -  ApBpSpJ^Qp  1 


I  32) 


where  vec(*}  denotes  vectorization  operator  (stacking  the  columns  of  a  matrix  on  top 

_  _i 

of  eacli  other),  respectively,  and  Qp  2  is  the  Hermitian  square  root  of  Q7 1 .  In  (1  :{:>), 
we  have  omitted  the  high-order  terms  of  {A,}  for  the  approximation. 

Hence,  for  a  large  number  of  data  samples,  the  optimization  problem  in  (  I  1 1) 
can  be  approximated  as: 


min  ||  vec(Qp  2  XII^)  -  vec{Q,,  A,,B,,S„)  ||J. 

B„ 

with  Bp  =  Diag(Bp|,  B„p).  I  13 


To  solve  the  above  optimization  problem,  we  introduce  below  two  partitioned 
matrix  operations  and  two  lemmas  without  proof  (see  [21]  for  the  detailed  proofs). 


lli 


Definition  1.  Let  G  be  a  partitioned  matrix  with  GltJ  being  the  (i,  j)f.h  (/,  / 

1,  2,  ■  •  •  J)  submatrix  of  G.  Then  the  block-diagonal  vectorization  operation  is  dr f ini  d 
by 

vecb{G)  =  [t/et;(Gi,i)r  uec(GL>ta)7'  •  •  •  vec(G,u)T)T .  i  I  3  1) 

Definition  2.  Let  £  and  fl  be  two  partitioned  maince s  with  conformal  pail/lmi  , 
and  with  £fiJ  and  J  being  the.  (i,  j)th  submatrices  of  £  and  S!.  r<  spn  tin  Ip  I'h,  t 
the.  generalized  Khatri-Rao  product  [i2j  is  defined  by 

[£  ©  H] irj  ^  £ij  <g>  fly,  (4-35) 

where  [■],.;  denotes  the  (i,j)th  submatrix  of  the  given  partitioned  matrix  and  @  denotes 
the  Kronecker  matrix  product. 

Note  that  both  the  block-diagonal  vectorization  vecb(-)  and  the  generalized 
Khatri-Rao  product  ©  are  defined  based  on  a  particular  matrix  partitioning,  i.e.. 
different  matrix  partitionings  will  lead  to  different  results.  Throughout  this  report 
the  partitioned  matrix  operation  are  all  based  on  the  partitionings  in  { 1  lb)  -  ( l  j  I  j. 
It  is  also  worth  pointing  out  that  matrix  partitioning  may  be  “inherited'’  through  nm 
trix  operations.  For  example,  for  the  partitioned  matrix  A  given  by  i  •  .  \!‘  A 
is  a.  p  x  p  partitioned  matrix  with  the  (i,  j)th  (i,j  —  1.  2  .  p\  siibm.tinx  ln'iir. 
Afl(§i)A(dj). 

Lemma  1.  Let  £  and  £2  be  two  partitioned  matrices  with  l\  bloc  I  rows  and  .1  him  I, 
columns,  and  let  G  be  a  block- diagonal  matrix  with  compatible  dimensions  and  con¬ 
formal  partitioning  with  £  and  SI.  Then 

vec6(£GnT)  =  (SI  ®  £)ner;6{G).  il  3b) 

Lemma  2.  Let  U  and  V  be  two  partitioned  matrices  with  1  block  row  ami  .1  blmt 
columns,  and  let  H  and  F  be  two  partitioned  matrices  with  1  block  row  ami  A  him  I 
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columns  and  with  compatible  dimensions  with  U  anil  V.  respectively.  Thin. 


(U®  V)"(H©F)  =  (U"H)@  (VWF) 


Now,  let 


i.e., 


(3p  =  vecb{Bp), 


vecr(Bf;[)  •■■vecT(B, 


-\T 


By  using  Lemma  1  (with  K  =  1  and  ./  =  p),  we  obtain: 

vector4  ApBpSp)  =  [Sj*  ®  {Q^A)]/3„  4  r,£, 

Hence, 

||  vec{Q^XnSjt)  -  vec(Qp4A„BpSp)  ||2 
=  ||vec{QJT'xnsj-r,A||2 


>  II  vec(Q;7^XIlgf))  ||2  -vecw(Q^Xnsfi)f,J(f'/fp)-|r(/,/v«  (QJ.  JXFI, 
-  L  tr( Q; 1 R)  -  LM  -  vech" ( A  " Q" 1 XS" )  x 


(M"f  ®  (Af  Q;1  Ap/T  vecb{A"Q-lXS^), 


where  we  have  used  Lemmas  1  and  2,  and  the  equality  holds  when 


0  =  (fpfpJ-'fpvectQ^XrijJ. 


By  using  (1  23),  (i  25),  (1  32)  and  (1  11).  it  follows  that 

1 


max/(X[Hp) 

Q.B,, 


(7T e)^g{0u  ,  Op)  |Q,_ 


(  I  :S7 1 

(4-38) 

(4  3‘)) 

i  1  in 

U 

(I  II) 
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where 


90i ,  =  1  —  M  +  tr(Q“‘R)  -  —  vecb" ( A" Q p  *XS" )  x 

[(SPS£')r  ®  (A^Q-'A,,)]"’  vecb(A"Qp'XS").  (4-44) 


Similarly,  we  have: 


max  /{X|Hp+i) 

.Q-Bp+I 


(7re)w  yiO.tii.  ■  ■  ■  .0,,}  q,,  , 

where  Q;,+i  arid  g(G,6 1,  ■  *  •  ,  0p)  are  defined  similarly  to  Qfi  in  [  '  ]  and  <y(H, 

in  (  I  I  I),  respectively. 

Substituting  ( I  hi)  and  (  I  15)  into  (  I  22)  yields  the  conditional  GLR: 


>(0|{%})  =  /]  ,  -0,,)  IQ,,+  1|  1 

1  Sf(ffi,”-,^)|Q„|  J 


{  I  Hi) 


4.3  Iterative  Generalized  Likelihood  Ratio  Test  (iGLRT) 

The  basic  idea  of  the  iterative  generalized  likelihood  ratio  test  (iGLRT)  is  to 
detect  and  localize  targets  sequentially.  In  each  step  of  the  iteration,  the  results 
from  the  previous  iterations  and  steps  are  exploited  for  the  detection  and  localization 
of  new  targets  by  calculating  cGLR.  Specifically,  we  first  perform  GLRT  to  get  the 
location  of  the  dominant  target,  and  the  following  targets  are  del  cried  and  h ><  ali/M 
by  using  cGLRT  conditioned  on  the  most  recently  available  estimates.  The  detailed 
steps  of  iGLRT  are  described  in  Table  1. 

Once  the  locations  of  the  targets  are  determined,  the  amplitudes  of  the  reflected 
signals  can  be  estimated  by  using  the  AML  estimator  in  [21]: 

0*  =  [(Agq*‘A*)  ®  (SgS^JT’wcbtAjQ^XSJ),  (4- 47) 

where  A^,  and  Q,--  are  defined  similar  to  Ap.  Ap  and  Q,,  in  (  I  is),  ( I  10)  and 
(1  27),  respectively. 


Table  -1  1 :  iGLRT  algorithm 


Step  I: 

-  Calculate  p(6)  in  (  I  12)  for  each  0. 

-  Compare  p(0)  to  a  threshold  pa.  if  p(0)  <  plt  for  all  II.  ihni  Sn>|. 
otherwise,  H\  =  arg max^ p{8).  go  to  Step  II. 

Step  II:  For  k  =  1,  2,  ■  •  ■ ,  do  the  following: 

-  Calculate  p{#|{^t}J=i)  in  ( 1  Hi)  for  each  0. 

-  If  <  pa  for  all  0,  then  go  to  Step  III; 

otherwise,  0b+[  =  arg  max#  p{0\{H,  ) . 

Step  III:  Repeat  the  following  substeps  until  convergence  for  k  —  1.  2.  .  h 

(suppose  that  K  targets  are  detected  in  Steps  I  and  II) 

*  Calculate  for  each  0. 

-  Update  Bk  by  argmaxp 


We  note  that  Step  III  of  the  above  iGLRT  algorithm  actually  mmiiuizi'>  i  h< 
function  g{6\,-  ■  ■  ,0%)  with  respect  to  by  using  the  cyclic  minimization  (CM  1 

technique  [33].  Under  a  mild  condition,  i.e.,  L  3>  NK.  we  have  <i(d\  ■  ■  ■  .  0h  t  n 
Furthermore,  we  know  that  the  CM  algorithm  monotonically  decreases  the  cost  func¬ 
tion.  Hence,  the  iGLRT  algorithm  is  convergent.  When  K  is  the  true  timrtbei  nl 
targets,  iGLRT  reduces  to  an  approximate  (parametric)  maximum  likelihood  rsnum- 
tor.  As  we  will  show  via  numerical  examples,  the  mean-squared -error  (MSI;  ul  i  In 
estimate  of  iGLRT  approaches  the  corresponding  Cramer-Rao  bound  (CRB)  fm  a 
large  number  of  data  samples.  On  the  other  hand,  we  note  that  iGLRT  needs  onl\ 
one-dimensional  search  and  hence  is  computationally  efficient. 


CHAPTER  5 

NUMERICAL  EXAMPLES 

In  this  chapter,  we  first  compare  the  Cramer- Rao  bounds  (CRBs)  for  MIMO 
radars  with  different,  antenna  configurations,  and  then  present  the  detection  and 
localization  performance  of  the  proposed  methods, 

5.1  Cramer- Rao  bound 

We  first  study  the  Cramer- Rao  bound  under  various  antenna  configurations 
Consider  a  MIMO  radar  system  with  M  =  N  =  8  antennas  for  transmitting  and 
receiving.  We  assume  that  the  receiving  and  transmitting  antennas  are  grouped  into 
multiple  subarrays  {each  being  a  uniform  linear  array  with  half- wavelength  spacing 
between  adjacent  elements): 

•  MIMO  Radar  A:  1  subarray  with  8  antennas  for  transmitting  and  receiving; 

•  MIMO  Radar  B:  2  subarrays  each  with  4  antennas  for  transmitting,  and  1 
subarray  with  8  antennas  for  receiving; 

•  MIMO  Radar  C:  8  subarrays  each  with  1  antenna  for  transmitting,  and  1  sub¬ 
array  with  8  antennas  for  receiving; 

•  MIMO  Radar  D:  2  subarrays  each  with  4  antennas  foi  transmitting  amt  m  m 
ing. 

We  assume  that  the  transmitted  waveforms  are  linearly  orthogonal  lu  each  uilin  .uni 
the  total  transmitted  power  is  fixed  to  be  1,  i.e.,  R„m,  =  ^ I 

We  consider  a  scenario  in  which  K  =  3  targets  are  located  at  B\  =  —  40°,  — 
—4°  and  0$  —  0°,  and  the  elements  of  {B^. }jj=1  are  independently  and  identicalh 
distributed  (i.i.d.)  circularly  symmetric  complex  Gaussian  random  variables  with 
zero  mean  and  unit  variance.  There  is  a  strong  jammer  at  It)  with  amplitude  Hill, 
i.e.,  40  dB  above  the  reflected  signals.  The  received  signal  has  L  —  128  snapshots 
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and  is  corrupted  by  a  zero- mean  spatially  colored  Gaussian  noise  with  an  unknown 
covariance  matrix.  The  (p,  r/)th  element  of  the  unknown  noise  covariance  matrix  is 

Figures  5  1(a)  and  5  1(b)  show  the  cumulative  densih  functions  (C 'Dl  -  ><|  i  h< 
CRBs  for  MIMO  radar  with  various  antenna  configurations  when  .SNR  20  illi  I  In 
CRB  of  62  is  similar  to  that  of  9 3  and  hence  is  not  shown.)  The  CDFs  me  obtained 
by  2000  Monte-Carlo  trials.  In  each  trial,  we  generate  the  elements  of  , 

randomly,  and  then  calculate  the  corresponding  CRBs  using  (A  IS)  given  by  in  the 
Appendix.  For  comparison  purposes,  we  also  provide  the  CDF  of  the  phased-array 
(single-input  multiple-output)  counterpart,  i,e.,  the  special  case  of  the  above  MIMO 
radar  when  N  =  1,  with  the  same  total  transmission  power.  As  expected,  the  MIMO 
radar  provides  much  better  performance  than  the  phased-array  counterpart  Due  to 
the  fading  effect  of  the  elements  of  {B0fc}|=!,  the  CRB  of  MIMO  Radai  A  vanes 
within  a  large  range.  Within  a  95%  confidence  interval  (i.e.,  when  CDF  varies  tmm 
2.5%  to  97.5%),  its  CRB  for  9]  varies  approximately  from  5  x  10~7  to  5  x  ]()  The 
CRBs  for  MIMO  Radar  C  varies  within  a  small  range. 

To  evaluate  the  CRB  performance,  we  define  an  outage  (  Tlli  lot  i 
probability  p,  denoted  by  CRB?J,  as: 

P(CRB  >  CRBp)  =  p.  1 5  1 1 

Figures  5  2(a)  -  (d)  show  the  outage  CRBo.m  and  CRB0  i  of  9X  and  9$,  as  functions  of 
SNR.  As  expected,  the  SNR  gains  depend  on  the  probability  p.  As  we  can  see,  when 
p  —  0.01,  MIMO  Radar  C  outperforms  the  other  radar  configurations,  and  provides 
around  20  dB  and  12  dB  improvements  in  SNR  compared  to  the  phased-array  and 
MIMO  Radar  A,  respectively.  On  the  other  hand,  Figure  5  2(d)  shows  that  MIMO 
radars  A  and  B  outperform  others  when  p  =  0.1. 
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5.2  Target  Detection  and  Localization 

We  focus  below  on  MIMO  Radar  B,  i.e.,  a  MIMO  radar  system  with  2  miImi 
rays  (each  with  4  antennas)  for  transmitting  and  1  subarray  (with  8  antennas)  lot 
receiving. 

We  first  consider  a  scenario  in  which  3  targets  are  located  at  9i  =  —40°,  9>  = 
—20°  and  9$  =  0°  with  the  corresponding  elements  in  BW| .  and  B#a  being  fixed  to 
2,  2  and  1,  respectively.  The  other  simulation  parameters  are  the  same  as  for  Figure 
5  1.  The  Frobenius  norm  of  the  spatial  spectral  estimates  of  B,,  versus  9,  ubr.iim-d  h\ 
using  LS,  Capon  and  APES  are  given  by  Figures  5  3{a)-(c).  For  comparison  purpo^o 
we  show  the  true  spatial  spectrum  via  dashed  lines  in  these  figures.  V-  sera  limn 
Figure  5  3.  the  LS  method  suffers  from  high-sidelobes  and  poor  resolution  pioliletus 
Due  to  the  presence  of  the  strong  jamming  signal,  the  LS  estimator  fails  to  work 
properly.  Capon  and  APES  possess  excellent  interference  and  jamming  suppression 
capabilities.  The  Capon  method  gives  very  narrow  peaks  around  the  target  lm  . < i  mi 
However,  the  Capon  estimates  of  B0,,  B02  and  Bg,  are  biased  down  waul,  [’be  \  I '  I  s 
method  gives  more  accurate  estimates  around  the  target  locations  hut  its  resolution 
is  worse  than  that  of  Capon.  Note  that  in  Figures  ’>  3(a)-(c),  a  false  peak  omits  at 
9  ~  10°  due  to  the  presence  of  the  strong  jammer.  Although  the  jammer  waveform  is 
statistically  independent  of  the  waveforms  transmitted  by  the  MIMO  radar,  a  false 
peak  still  exists  since  the  jammer  is  40  clB  stronger  than  the  weakest  target  and  the 
number  of  data  samples  is  finite.  Figures  5  3(d)- (e)  give'  the  CLHT.  and  the  it  >1.1;  I 
results,  as  functions  of  the  target  location  parameter  9.  For  convenience  in  limn > 
5  3(e),  we  have  included  all  cGLR  functions  obtained  by  iGLHT.  each  mdi<  atmv  .  >tn 
target.  As  expected,  we  get  high  GLRs  (cGLRs)  at.  the  target  locations  ami  low 
GLRs  (cGLRs)  at  other  locations  including  the  jammer  location.  By  comparing  t  he 
GLR  with  a  threshold,  the  false  peak  due  to  the  strong  jammer  can  be  detected  and 


rejected,  and  a  correct  estimate  of  the  number  of  the  targets  can  be  obtained  by  bot  li 
methods. 

Next,  we  consider  a  more  challenging  example,  where  82  is  —4°  while  all  the  other 
simulation  parameters  are  the  same  as  before.  As  shown  in  Figure  '<  ,((■).  1  lie  APES 
Capon  and  GLRT  methods  fail  to  resolve  the  two  closely  spaced  taig',i-  at  <> 
and  83  =  0°.  On  the  other  hand,  iGLRT  gives  weil-resolved  peaks  at,  mm  I  1  In  nn< 
target  locations,  To  illustrate  the  procedure  of  the  iGLRT  algorithm  wc  give  tin 
GLR.  and  cGLRs  obtained  in  Steps  I  and  II  of  iGLRT  in  Figures  .(a)-(dj  l-'iguies 
■r>  5(a)-(b)  show  tlie  GLR  p(9)  and  the  cGLR  p{8\8\).  respectively,  where  9\  is  the 
estimated  location  of  target  1  from  p{9).  As  we  can  see,  there  is  no  peak  at  around 
=  0°  in  both  figures.  Yet  a  clear  peak  is  shown  iri  p{8\8\.8-i)  in  Figure  "  'i(c). 
which  indicates  the  existence  and  location  of  target  3.  The  cGLR  p{8\8\ .  8>.  9:i)  in 
Figure  5  5(d)  shows  that,  110  additional  target  exists  other  than  the  targets  at  8 8> 
and  83.  In  other  words,  the  iGLRT  method  correctly  estimates  the  number  of  targets 
to  be  3, 

Now  we  consider  the  elements  in  B0],  Bg2  and  Bw,,  as  i.i.d  complex  Gaussian 
random  variables  with  mean  zero  and  unit  variance.  The  other  parameters  are  the 
same  as  those  in  Figure  5  (>,  The  Figures  5  tj(a)--(b)  present  the  (’l)Fs  nl  ilu-  MSI 
of  8\  and  83  as  well  as  the  CRBs,  when  SNR  =  20  dB  and  L  —  128.  As  we  can 
see,  the  MSEs  of  the  iGLRT  are  very  close  to  the  corresponding  CRBs.  Figures 
5  7(a)-(b)  show  the  outage  MSE0.i  and  CRBfu  when  p  -  0,1  as  functions  of  SNR 
when  L  —  128.  Again,  the  MSEs  me  very  close  to  the  corresponding  the  CRBs. 
and  decreases  almost  linearly  as  SNR  increases.  Figure  5  g  gives  the  outage  MSE0.1 
and  CRBn.i  as  functions  of  L  when  SNR=20  dB.  As  expected,  the  outage  MSEn  1 
approaches  the  corresponding  CRBn.i  as  L  increases. 
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(a)  0i  (lij  0t 

Figure  5  1:  Cumulative  density  functions  of  the  Cramer- Rao  hounds  for  (at  0\  mu  I 

(b)  e-s, 


(a)  CRBo.m  for  0t 


(1))  CRBo.m  for  02 


Figure  5  2:  Outage  CRB  versus  SNR.  (a)  CRBo.m  for  l> |.  (h)  CRBum 
CRBn  i  for  0],  and  (d)  CRBo.i  for  ()■>. 
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(e)  iGLRT 


(d)  GLRT 


Figure  5-3:  Spatial  spectra,  and  GLR  and  <:GLR  Pseudo- S | n i  in  .when  D  n 

02  =  -20°,  and  &3  —  IF.  (a)  LS.  (b)  Capon,  (c)  APES,  (d)  GLUT,  mid  u  i(  d.K  I 
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Figure  5-4:  Spatial  spectra,  and  GLR  and  cGLR  Pseudo-Sport  ra.  when  0,  lit 
9%  =  -4°,  and  0;t  =  0°.  (a)  Capon,  (b)  APES,  (c)  GLRT.  and  (.1)  iCLlH 


(c)  p(6\6\,Qt)  (cl) 


Figure  5-5:  GLR  and  cGLR  Pseudo-Spectra  obtained  in  Steps  l  and  II  of  iCLRT. 
when  0 ,  = -40°,  02  =  -4°,  and  0:I  =  0*.  (a)  p{9)t  (b)  p(6\ffi),  (c)  p(Wh.fh).  mi  (d) 
p{Q\Q\ .  0i,  0:\ ) . 


Figure  5-6:  Cumulative  density  functions  of  the  CRBs  and  MSEs  lui  mi  n,  ,ind  t 

h- 


MS£ 
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Figure  5  7:  Outage  CRBn.j  and  MSEn.j  versus  SNR  for  (a)  and  (b)  0:x. 
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(a)  9X  (It)  <h 


Figure  5-8:  Outage  CRBo  i  and  MSEo  i  versus  number  of  data  samples  tni  im  w,  .md 
(f>)  03- 


CHAPTER  6 
CONCLUSIONS 

We  have  considered  a  multiple-input  multiple-output  (MIMO)  radar  system  with 
a  general  antenna  configuration  that  can  be  used  to  achieve  both  the  coherent  pro¬ 
cessing  gain  and  the  spatial  diversity  gain.  We  have  first  introduced  several  spatial 
spectral  estimators,  including  Capon  and  APES,  for  target  detection  and  paiannin 
estimation.  By  using  our  results  on  the  growth  curve  models,  we  have  provided  a 
generalized  likelihood  ratio  test  (GLRT)  and  a  conditional  generalized  likelihood  ra¬ 
tio  test  (cGLRT),  and  then  proposed  an  iterative  GLRT  (iGLRT)  procedure  lor  the 
MIMO  radar  system.  Via  several  numerical  examples,  we  have  shown  t  hat  t  In-  it  il.li'l 
method  can  provide  excellent  target,  detection  and  parameter  estimation  performance 
at  a  low  computational  cost. 
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APPENDIX  A 
CRAMER.- RAO  BOUND 

Consider  a  MIMO  radar  system  with  K  targets.  Then  the  received  smn.d  >  .n. 
be  written  as 

K 

X  =  53  M0k)Bll,VT{0k)$  +  Z  (A  I) 

A  =  i 

Let 

0  =  [Sj  •  •  •  tin] 1 ,  i  A  1 1 

0  =  [vecr(Be,)  •  ••  vec7’(BjK)]7'(  (A  3) 

and 

0n  =  Re{0)  (3,  =  Im(/3),  (A-4) 


where  Re(')  and  Im(-)  denote  the  real  and  imaginary  parts,  respectively.  Assume  that 
the  columns  of  Z  are  i.i.d.  circularly  symmetric  complex  Gaussian  mm  him  ve<  i<u> 
with  zero-mean  and  an  unknown  covariance  matrix  Q. 

Using  the  same  argument  as  in  Appendix  A  of  .  we  know  tlmi  t  he  unknown.- 
in  Q  will  not  affect  the  CRBs  of  0  and  0.  Hence,  we  need  only  to  calculate  the 
following  Fisher  information  matrix  with  respect  to  0,  0n  and  0,\ 


PIM 


F(0.0)  F  (0,0R)  F(0.j3,) 

F (0R,0)  F(/3/f, 0R)  F(0r,  0t) 
F  (0„8)  F(0h0R)  F(/3„/3 ,) 


(A  5) 


where  F(a,/3)  denotes  the  Fisher  information  matrix  with  respect,  to  a  and  0, 
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Note  that 


F(8„  dj)  =  2  Retr  < 


ZLM9k)B0kVT(dk)$ 


i  a 


08, 


Q 


08, 


(A  (i) 


and 


£LA(4)B,4v^.)$ 


08, 


1  =  M&i)B0kyT(8t)^  +  A(8,)Bt)y  (8,  )<h  (A  7. 


where 


AW  =  ^i  and 


08 


08 


(A  rt) 


Inserting  (A  7)  into  (A  (i)  and  after  some  matrix  manipulations.  we  obtain. 

F{0i,8j)  =  2L  Re  tr  j  [A 11  (8, )  Q  1  A(%}]  ;B„  V'(W,)R,MV(«  lU  f]} 

2LRetr{[A"{01)Q"1A(^)](B^V7'(^)R,l,<l,V*(^)B'']}  + 

J  (A  'M 

2LR(!tr{[A,/(^;)Q-1A(0J)]jBw  VT(^)R(N,V*(y,)B;;]}  + 
2LRetr{[Aw(0i)Q-,A{%|[B%Vr(^)R*<l.V*(4)B^]}, 
where  R,m>  =  is  the  covariance  matrix  of  the  transmitted  waveforms.  Hence, 


F{0.  9)  =  2  R.e(Fo</),  (A  10) 

where  Fo»  is  a  K  x  K  matrix  with  its  (i,  j)  element  being: 

[Fesh  =  Ltr{[A(^)y'Q_1A(%)]  [B0] VT(8J)Rt] . V‘(0,)B"]}  + 

Ltr  { [A"(0,)Q-' Affl,-)]  [B0V7'(^)R<N,VT{0,)B//]}  + 

J  (A  11) 

L  tr  |  [AW(6),:)Q-1  A(0j)]  [B0  (Vr(0J)RM,V-(fil,  )B,"] }  + 
LRe{[A,/(ei)Q-lA(^)][B0,VT{6)J}R+4.V*(#()B^j} 


Similarly,  we  have 


FT(0,/3/{)  =  F(/3li.fl)=2Re(Fao)t  (A  12) 

FT(0./3,)  =  F(/3,,0)  =  2Im(Fw).  i  \  !  (. 


—  F(/3/,/3/)  —  2 Re(Ffltf),  \  111 

and 

F  Tt0Rtpr)  =  F((3h0R)  =  2Inx(Fw),  (A  15) 

where  F0/J  and  F$$  are  both  partitioned  matrices  with  K  x  K  and  K  x  K  blocks  and 
with  their  submatrices  being,  respectively. 


[F»]«  =Lvec{iA"(0,)Q-1A(^)]  [B0j VT(0J)R4><I,V(ffr)j}  + 

Lvec  j[Aw(0i)Q’1A(0j)]  [BflVr(0J)ILM,V*(01)]} 


(A  Hi) 


and 


My  =  L[VT(ft)R.M>V*(^}]  ®  [A^fftJQ-1  A(0,)],  (A  17) 


with  ®  denoting  tiie  Kronecker  product. 

Substituting  Equations  (A  !))  -  {A  15)  into  (A  U),  and  alter  some  nmlri.s  m.i 
nipulations,  we  get 


CRB(ff)  =  l-{  Re  [Foo  -  FWF^F"]  } 


(A  IS) 


and 


CRB(£»)  =  F^  +  F^F^CRB(O)  FMF;‘. 


I A  Id  i 
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